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By means of atomistic molecular dynamics simulations, we study cholesterol-DPPC (dipalmitoyl 
phosphatidylcholine) bilayers of different composition, from pure DPPC bilayers to a 1:1 mixture 
of DPPC and cholesterol. The lateral-pressure profiles through the bilayers are computed and 
separated into contributions from the different components. We find that the pressure inside the 
bilayer changes qualitatively for cholesterol concentrations of about 20% or higher. The pressure 
profile then turns from a rather flat shape into an alternating sequence of regions with large positive 
and negative lateral pressure. The changes in the lateral-pressure profile are so characteristic that 
specific interaction between cholesterol and molecules such as membrane proteins mediated solely 
via the lateral-pressure profile might become possible. 



I. INTRODUCTION 

From a macroscopic point of view a planar lipid mem- 
brane in equilibrium has, by definition, a vanishing sur- 
face tension. However, on a microscopic level there is, 
within the membrane, local lateral pressure, i.e., pres- 
sure tangential to the interface. Only when all of the lo- 
cal contributions are summed and averaged along the bi- 
layer normal, the net pressure vanishes. The local lateral 
pressure arises from the different structural components 
of the lipid bilayer l|MarsrA Il99fil) and, almost counter- 
intuitively, each of these contributions can be of the order 
of several hundreds of bars. 

The importance of lateral pressure profiles ha s been 
discussed in several recent reviews llBen-Shauli Il995t 
| BezrukovLl20MlEckenholll200ltl^nnunenLl200(ilMarshl 
1996). The increased interest in understanding lat- 
eral pressure profiles is easily understood as lateral 
pressure has been proposed to have an import a nt role 
in, for example , general anaesthesia ijCantoit Il997al 
lEckenho"fiL l200l |) or inh i bition and regulation of pro- 
tein function llde Kruiiffl. Il99 7|) (see referenc es 25-50 in 
Ref . Ivan den Brink- van der Laan et all 120041 for a list of 
proteins for which a relation to lateral pressure has been 
suggested). 

Due to the absence of good probes for lateral pres- 
sure, direct experimental measurements are difficult 
and onl y a single experime ntal study exists at the 
present ijTempler et all Il998|) . In contrast, computer 
simulations have for the past ten years been able 
to supply direct, yet not straightforward, access to 
study pressure profiles and their response to changes 
in the membrane. Calculations of lateral pressure 
profiles are a delicate matter as already the intro- 
duction of small simplifications to the system can 
render the results questionable. For example, the 
reported results for pure lipi d bilayer systems from 
coarse-grained simulati o ns IIGoetz an d LiDOwskv, 1998; 
Harries and Ben-Shaul Il997t IShillcock and Lipowskv. 
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In this paper, we study the effects of cholesterol 
on lateral pressure. Cholesterol is an essential com- 
ponent of all Eukaryotic cell membranes where it 
plays a crucial role for both stat i c structure an d dy- 
namics ( bimons and Ikonenl l2000t lYeaglet Il985|) . In 
partic ular it regulates the fluid i ty of the cell mem- 
brane ijMcMullen and McElhanevL ll996). All of this goes 
along with changes in the lateral pressure profile. Indi- 
rect evidence for the importance of cholesterol on the lat- 
eral pressure comes from studies of membrane channels 
which are hig hly sensitive to th e pressure of their en- 
viron ment ijHamill and Martinad l200ll ISukharev et all 
119971) . One such channel is the nicotinic acetylcholine re- 
cepto r which ceases to fun ction in the absence of choles- 
terol ijRankin et allll997|) . 

Apart from the fundamental importance of cholesterol 
and its effect on lateral pressure, the study of choles- 
terol is interesting for a second reason: cholesterol is a 
highly specific molecule. Already small modifications of 
its sterol structure l ead to significa n t changes of the mem- 
brane properties ijEndress et all l2002t IScheidt et all 
2003). Generic theories of lateral pressure describe a 
molecule basically only by its volume and its rigidity. 
When such theories are applied to cholesterol l)Cantorl 
1999b), they thus inevitably fail to capture many es- 
sential features of cholesterol, and even predict pressure 
changes of wrong sign in some parts of the bilayer. 

In this article, we thus use atomistic molecular dy- 
namics simulations to study lateral pressure profiles for 
six different systems, ranging from a pure DPPC bilayer 
to a bilayer consisting of a 1:1 mixture of DPPC and 
cholesterol. This paper is, to the author's knowledge, 
the first detailed atomistic computational study address- 
ing the effect of cholesterol or other small molecules on 
the build-up of the lateral pressure profile in phospholipid 
membranes. 
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II. LOCAL PRESSURE 

The pressure tensor p can be computed as 

P=2E-V, (1) 

from the kinetic energy density tensor E and the config- 
uration stress tensor E. Both of the latter quantities can 
be expressed in terms of atomistic positions, velocities 
and forces as 

E = - ^2 m iVi <8> Vi , (2) 
i 

E = 1 4- ® r y - , (3) 

and are thus accessible in a MD simulation. While 
the above expressions are, strictly speaking, defined 
only if the summations are extended ov er the en- 
tire simulation volume l)Heinz et all [2004) , in prac- 
tise it is possible to divide the different contributions 
into sl ices according to the posit ions of the involved 
atoms ijLindahl and EdhohnLl2000|) . 

This implies, however, that the force Fij between 
particles i and j is known explicitly. This is not the 
case if a multipole or lattice based method (such as 
PME) is used to evaluate electrostatic interaction. In 
the first reported atomistic compu tation of pressure 
profiles (|Lindahl and Edholml [2000) electrostatics was 
therefore truncated a distance of 1.8 nm. As it is known 
by now, however, using abrupt truncation, especially at 
such a short dist ance, introduces sign i ficant artifacts into 
bilayer systems l|Anezo et all l2003t iPatra et all 120031 
l2004albf) and thus needs to be avoided. We thus use 
reaction-field technique that has been shown to give re- 
sults consi stent with the appli cation of long-range elec- 
trostatics l|Patra et al 1 1200411 while at the same time 
employing explicit expressions for . 

The global pressure is equal to the average of the lo- 
cal pressures, and any condition on the global pressure 
thus translates onto the pressure profile. Since the out- 
side of a bilayer is at equilibrium with the environment 
(i.e., approximately 1 bar in most cases), the average lo- 
cal pressure has to be equal to that value. If the global 
pressure would be different, the system would react by 
shrinking or expanding, and thus would not be in equi- 
librium. 

Still, there is nonvanishing local pressure even in equi- 
librium. The existence of an interface between the water 
and the lipid goes along with an energy penalty which 
could be lowered by packing the bilayer more densely, 
thereby decreasing the area per lipid. Steric constraints 
between the lipid tails prevent this from happening. The 
equilibrium value of the area per lipid is thus a com- 
promise between the "wishes" of the head groups and 
the tails. This is directly reflected in the lateral-pressure 
profile. Since the interface region prefers a further reduc- 
tion of the area per lipid, the local lateral pressure there 



is negative (pointing inwards) whereas it is positive in 
the tail region (pointing outwards). 

The lateral pressure profile in equilibrium thus is a di- 
rect reflection of the inhomogeneity of the bilayer along 
the bilayer normal. In contrast, the bilayer is homoge- 
neous parallel to the bilayer interface. In equilibrium, 
the normal component of the local pressure thus has to 
vanish everywhere. (More correctly, it has to be con- 
stant and equal to the applied external pressure of 1 bar.) 
In contrast to the above theoretical argument, a non- 
vanishing normal pressure component is found in the nu- 
merical simulations, in magnitude about 10 % of the lat- 
eral component. We attribute this to the use of distance- 
constraints in our simulation. For computational effi- 
ciency, the distance between bonded atoms is kept con- 
stant, meaning that the system cannot locally expand or 
contract to reduce the loca l pressure. This "con s traint 
force" was studied in Ref. iLindahl and Edholm, 2000, 
and the values quoted there are able to explain the val- 
ues of the normal pressure found in our simulation. Since 
effects of numerics are isotropic on average, we can thus 
improve the lateral component of the pressure by sub- 
tracting the normal component from it before the analy- 
sis. This will be done throughout this paper. 



III. SIMULATION DETAILS 

We study lipid bilayers comprised of 128 molecules 
(64 per leaflet), at various ratios of dipalmitoyl phos- 
phatidylcholine (DPPC) and cholesterol, hydrated by 
3655 water molecules. D PPC molecules are de- 
scrib ed by the model from Ref. iTieleman and Berendserl 
Il996t which util i ses t he description of lipids from 
Ref. iBereer et ail 1199 7 } Choleste r ol wa s described by 
the m odel from Ref. iHoltie et all l200ll and the SPC 
model (jBerendsen et all Il98lj) was used to describe wa- 
ter. The simulations were performed usin g the Gromacs 
packa ge, both in the standard release ifLindahl et all 
1200 ID and in an adapted version that allows the com - 
putation of local pressures ijLindahl and Edholml |2000|) . 

For computing el ectrostatics interacti ons, we employed 
a twin-range setup l)Bishop et allll997|) in which the in- 
teractions within a distance rn st — 1.0 nm were eval- 
uated at every integration step, and those between rn st 
and r cu t = 2.0 nm only eve ry tenth integration time step. 
A reaction-field approach ijTironi et all Il995|) was used 
to account for interaction outside of r cu t by assuming 
a homogeneous dielectric with e — 80. Lennard- Jones 
interaction was truncated at 1.0 nm. 

DPPC, cholesterol and water molecules were sepa- 
rately coupled to a heat bath at temperature T = 323 K, 
and the pressure was kept at 1 bar both using the 
Berendsen algorithms l)Berendsen et allll984j) . The size 
of the simulation box in the plane of the bilayer (a;- 
y plane) was allowed to fluctuate independently of its 
height. 

As initial configurations for all simulations we used the 
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final c onfigurations of 100 ns simulations ijFalck et all 
2004b) where el ectrostatics were treated by particle-mesh 
Ewal d (PME) i|Essmann et all 119951 iFrenkel and Smhi 
|2002|) . The actual simulation procedure was divided into 
two steps. In a first step, 50 ns trajectories were gener- 
ated from these structures for all cholesterol concentra- 
tions using the standard version of Gromacs . Electrostat- 
ics w as treated by reaction field technique ijTironi et all 
1995), and the generated configurations were saved ev- 
ery 10 ps. The bond lengths of DPP C molecules were 
constrained by the LINCS algorithm ijHess et all Il997|) 
and wate r molecules were kept rig i d by t he SETTLE al- 
gorithm ijMivamoto and Kollmanl E?92), such that an 
integrator time step of 2.0 fs could be used. In the 
second step, the pressure profiles were generated from 
this trajectory using an adapted version of Gromacs (see 
below for details). In the se runs, the SHAKE algo- 
rithm (|Rvckaert et allll977|) was used to constrain bond 
lengths 

Resolving the pressure spatially makes the simulations 
about one order of magnitude slower than a normal sim- 
ulation. The pressure profiles were thus computed in a 
second step, based on the saved configurations. Starting 
at every saved configuration, a 4 ps simulation was run to 
compute the pressure profile. The first 2 ps were ignored 
to rule out effects of finite precision of the saved configu- 
rations. Unless mentioned otherwise, only the final 30 ns 
of each simulation were included in the analysis. 

The instantaneous pressure p, computed from Eq. 
is fluctuating quickly in time. Even when the instanta- 
neous spatial average over the entire simulation box is 
considered, the pressure easily changes by several hun- 
dred bars within a single integration time step. Com- 
puting a statistically relevant pressure profiles thus is 
numerically challenging since one has to sample a large 
number of configurations. We evaluated the pressure pro- 
file for a number of simulation frames that is far larger 
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FIG. 1 Temporal development of the area per lipid for the 
different cholesterol concentrations from t = ns to t = 50 ns. 
The results computed using PME are shown in condensed 
form from t = —100 ns to t — ns. 



than in previous studi es where pressure profiles of bi- 
layers were computed JGulling-srud and SchultenL l2004t 
iLindahl and Edhoiml. 12000). This allowed us to divide 
the simulation box into 150 bins for computing the pres- 
sure profile, with the result having only negligible nu- 
merical noise. The remaining uncertainty in the pressure 
profile is mainly due to the temporal change of the area 
per lipid. To arrive at this data quality, a total of ap- 
proximately 25 000 hours of cpu time was needed. 

In previous MD studies on cholesterol-DPPC bi- 
layers, electrost a tics w ere handle d either by plain 
cutoff llTu et all Il998l) or PME ijFalck et all l2004bt 
iHofsafi et all l2003). In our simulations we used reaction- 
field technique, motivated by the mutual consistency be- 
tween reaction-fi eld technique and P ME found for pure 
DPPC systems l|Patra et all l2004bh . From our sim- 
ulations we found that the results with reaction-field 
technique and PME are almost identical also for mixed 
cholesterol-DPPC bilayers. We thus refrain from repro- 
ducing the entire standard set of quantities that are used 
to characterise a bilayer. Rather, we show only the tem- 
poral development of the area per lipid in Fig. HI T he re- 
sults of the earlier simulations ijFalck et all l2004b() done 
using PME are also shown there in compressed form from 
t = —100 ns to t = ns. From Fig.^it is then immedi- 
ately obvious that changing the electrostatics treatment 
from PME to reaction field at t — ns has no relevant 
effect onto the systems, and that the systems are in equi- 
librium. 



IV. PRESSURE PROFILES 

The computed lateral pressure profiles are summarised 
in Fig. [3 We first want to discuss a few general features. 

First, for low cholesterol concentration, the pressure 
does not decrease to zero at the edges of the simula- 
tion box. This is a sign that the bilayer is not fully 
hydrated, and it was estimated that an additional 4— 
5 water molecul es per lipid would be need ed for com- 
plete hydration ijLindahl and Edhoiml [2000]) . This num- 
ber is in agreement with our results. Increasing choles- 
terol concentration means a decrease of the number of 
lipid molecules, hence an increase in the number of wa- 
ter molecules per lipid. At 12.5 % cholesterol, the number 
of water molecules per lipid has increased by 4. For this 
and higher cholesterol concentrations the pressure in the 
bulk water phase indeed becomes zero, as can be seen 
from Fig. H 

Second, the lateral pressure profile is not flat but some 
parts of the bilayer would like to expand, at the same time 
that other parts would like to contract. While the net 
pressure, averaged over the entire bilayer, is small, the 
local pressure can be much higher than typical macro- 
scopic pressures. A simple estimate shows that the lo- 
cal lateral pressure in the bilayer core can be well over 
300 bar and cl ose to the interface can reach valu es exceed- 
ing 1000 bar ijGullingsrud and Schultenl l2004|) . Our re- 
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FIG. 2 Lateral pressure profiles through the bilayer for dif- 
ferent cholesterol concentration. The different curves have 
been shifted for clarity. The labels on the graphs give the 
cholesterol concentration. 



suits for pure DPPC agree b oth qualitatively and quanti- 
tatively with earlier results (Lindahl and Edholm, 2000) 
with small differences in the peak positions close the in- 
terface due to the different treatment of electrostatics in 
that study. 

Third, the magnitude of the local lateral pressure be- 
comes higher as the cholesterol concentration increases. 
While for small cholesterol concentrations the local pres- 
sure is of the order of few hundred bars, it increases to 
thousands of bar for the highest cholesterol concentra- 
tions. There is no straightforward explanation for this 
but this phenomenon is very likely related to higher bi- 
layer rigidity at high cholesterol concentration since pres- 



sure gradients are ultimately related to the elastic mod- 
ulus. 

Finally, the pressure profiles possess additional struc- 
ture in the presence of cholesterol. Without cholesterol, 
the lateral pressure in the lipid tail region of the bilayer is 
strictly nonnegative. Already for 5 % cholesterol, a small 
region of negative lateral pressure is seen in the figure. 
For cholesterol concentrations of 20 % and higher, addi- 
tional structure is seen. This reflects that cholesterol is 
not some generic structureless object but rather possesses 
an internal structure. One could call this the specific ef- 
fect of cholesterol, compared to the unspecific effects that 
are also observed. 



V. PARTIAL PRESSURES 

The pressure tensor Q arises from inter- and in- 
tramolecular pairwise forces [cf. E from Eq. (J3J] as 

well as from the kinetic motion of the atoms [cf. E 
from Eq. J2J)]. There are three components in the bi- 
layer system (DPPC, cholesterol and water), which gives 
six possible combinations for the pairwise forces, and E 
can thus be split into six contributions (DPPC-DPPC, 
DPPC-cholesterol, and so on), according to which kind 
of molecules the two atoms i and j causing the force Fij 
belong to. 

Similarly, the kinetic energy tensor E can be split into 
three contributions. Thus, we divide the lateral pressure 
profile p{z) into nine contributions - six from E and three 

from E. The partial pressure profiles in Fig. [3] show the 
different contributions marked by different colours. The 
size of the coloured area directly gives the pressure due to 
that contribution. The arrangement of the colours, i.e., 
whether a given colour is close to the z-axis or further 
away, has no physical meaning. The kinetic contribu- 
tion from the DPPC molecules is marked with the same 
colour tone but somewhat lighter (i.e., light red) than 
the DPPC-DPPC contribution from the forces (i. e., dark 
red), and so on. 

There exists a fundamental difference between the ten- 
sor E, in that it depends only on positions and forces and 
thus describes only the static properties of the system, 
and E, which depends only on the velocities and thus is 
a purely dynamic quantity. The kinetic energy per par- 
ticle in Eq. @ can be larger than k-gT as, e.g., bond 
vibrations or nonbonded interactions contribute energy 
to the translational degrees of freedom. The increase 
in kinetic pressure E depends on whether these interac- 
tions are mainly oriented parallel to the x-y plane or nor- 
mal to it. It thus complements the traditional deuterium 
or NMR order parameter |Scd| that is used to quantify 
orientational orientation along a chain ijTieleman et alj , 

E223). 

Each contribution can be either expanding (positive 
pressure) or contracting (negative pressure). The sum 
of all expanding contributions is shown as black line in 
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FIG. 3 Partial pressure profile through the bilayers. Expanding (p > 0) and contracting (p < 0) contributions are shown 
separately. The colours code the origin of the pressure; lighter colours stand for kinetic contributions. The label on the curves 
gives the cholesterol concentration (bottom). On the top, an enlarged version of the profile for 30% cholesterol is shown where 
important peaks are labelled according to which position in a DPPC and / or cholesterol molecule they relate to. Schematic 
drawings of these molecules are included on the right. The peaks for other cholesterol concentrations could not be labelled for 
space reasons but they are easily mapped to the top figure. 



the figure. The same applies to the sum of all contract- 
ing contributions. The difference between these two val- 
ues gives the pressure profile that was shown in Fig. [21 
It should be noted that this difference is significantly 
smaller in magnitude than each of the two original terms. 

As more cholesterol is added, the positions of the pres- 
sure profile peaks shift. Simultaneous analysis of both the 
pressure profile and the atom density profiles as a func- 
tion of cholesterol concentration allows us to find corre- 
lations between pressure and atom positions. The result 
of this analysis is included in Fig. [31 Labels "A"-"E" 
refer to peak positions caused by cholesterol, and labels 
"W"-"Z" to DPPC. 

Let us now take a closer look at the different regions 
in the bilayer. 



A. Interface region 

There is hardly any contact between water and choles- 
terol molecules, so one might be tempted to assume that 
the interface region (also referred to as headgroup region) 
of the bilayer should be hardly influenced by cholesterol 
at all. The opposite is the case, as is seen from Fig. [31 
Cholesterol reduces the area per lipid of the bilayer, and 
that strongly influences also the interface region, but the 
effects there are mainly generic and unspecific - any sub- 
stance that would reduce the area per lipid would have a 
similar effect. 

The outermost part of the bilayer is formed by the 
choline group of the DPPC molecules (labelled as "W"), 
and the total pressure there is negative. The attrac- 
tion is dominated by the solvation energy between water 
and the polar lipid headgroups. Most of the the con- 
tracting pressure is compensated by the positive pres- 



6 



sure among the water molecules. Entropy plays an im- 
portant role in the latter as water molecules become or- 
dered in the electrostatic field f rom t he zwitterionic head- 
groups (Linda hl and Edholml 12000). The behaviour of 
this region of the pressure profile depends only weakly 
on cholesterol concentration as the hydra tion energy de - 
pends only weakly on the area per lipid ijMarshl . Il996fl . 
Note that both the DPPC bilayer (red) and the water 
phase (blue) would like to expand, and it is only the inter- 
action between these two (green), commonly called sim- 
ply "hydrophobicity" , that is trying to reduce the area 
of the interface (green). 

A bit further down in the bilayer, the phosphate groups 
of the DPPC molecules are located. The mutual ar- 
rangement between the choline and p hosphate groups de- 
pends heavily on t he ar ea per lipid ijFalck et al.i r2004b; 
iGurtovenko et all l2004|) . For large area per lipid, and 
thus for low cholesterol concentration, both groups are, 
on average, located in the same plane. In this case, the 
effects of the phosphate and the choline group cannot 
be separated in the pressure profile. At sufficiently large 
cholesterol concentration, the choline group is tilted up, 
and an additional peak ("X") thus appears. The pres- 
sure is positive as the water molecules are oriented ac- 
cording to the charge of the choline group and thus have 
an unfavourable interaction with the oppositely charged 
phosphate group. 

The inner boundary of the interface region is marked 
by the sn— 3 carbon (labelled as "Y"). This is the furthest 
extension to where there is noteworthy penetration of 
water into the bilayer. Furthermore, this is the outermost 
position of cholesterol, occupied by its hydroxyl group 
("A"). Only inside this peak, cholesterol has, in addition 
to its change of the area per lipid, also direct effects. 



B. Acyl chain region 

The acyl chain region is of special importance to the 
study of cholesterol-containing bilayer as here not only 
the acyl chains of DPPC are located (hence the name) 
but also the rigid four-ring structure of the cholesterol 
body. In the presence of cholesterol, a characteristic 
structure in the pressure profile develops, to the point 
that one might speculate on whether this is sufficient to 
explain specific effects, such as observed for the nicotinic 
acetylcholine receptor l)Rankin et all Il997|) . 

Upon increase of cholesterol concentration, an order- 
ing of the DPPC acyl chains takes place. This ordering 
increases the attraction among the DPPC tails, thereby 
increasing the compressing pressure component, while at 
the same time introducing orientational correlations into 
the motion of the atoms of the tails, thereby increas- 
ing the expanding kinetic pressure. The kinetic pressure 
per DPPC molecule increases monotonously by almost a 
factor of 4 from % to 50 % cholesterol, and per choles- 
terol molecule it increases by a factor of 3. Cholesterol 
thus induces significant correlations into the motion of 



the atoms, and the ordering induced by cholesterol is 
reflected much stronger in the correlated motion of the 
atoms than it is in the order parameter |Scd|- 

The rigid four-ring structure of cholesterol of- 
fers a much more interesting chemical structure 
than the basically linear chains of the DPPC 
molecules. This shows also in the pressure profile 
that exhibits several well-pronounced peaks, mak- 
ing the pressure profiles of cholesterol-containing 
membranes very di fferent from previously re ported 
pressure profiles llGullingsrud and SchultenL l2004t 
iLindahl and Edholm! 120001) . There are several atoms in 
the cholesterol molecules that are correlated with peaks 
in the pressure profile, see the labels in Fig. |3| Steric 
arguments can explain only same of the correlations. 
This should come as no surprise since the pressure p(z) 
at given depth z in the bilayer cannot be computed 
solely from atomic information from that depth - the 
stresses in different parts of the membrane are coupled 
due to finite elastic modulus. The picture becomes even 
more difficult if one accounts also for chemistry, i.e., 
for favourable and unfavourable interactions between 
certain pairs of atoms. 

The orientational correlations in the kinetic pressure, 
related to entropy and giving rise to kinetic pressure, 
complicate matters even more. There is thus little point 
in trying to explain the origin of every peak in the lateral 
partial pressure profile, especially since at the moment 
there is little possibility to prove or disprove any sugges- 
tion. What is sure, however, is that all the features of 
the lateral pressure profile arc in some way related to the 
chemistry and structure of cholesterol. 

Cholesterol is by far the most prominent steroid in any 
Eukaryote but it is only one member of a large group 
of chemically related substances that are varying only 
by seemingly minor changes in their structure, such as 
a single additional double bound or a different place- 
ment of the two carbon atoms that are sticking out from 
the rigid cholesterol body. Still, their pr operties in the 
membrane are distingu ishably different ijEndress et all 
l2002t [Scheidt et all 12003^ . and only in special cases one 
type of steroid can replace another in biological systems. 
Our results suggest that cholesterol relatives such as, 
e. g., lanosterol or ergosterol, would have a pressure pro- 
file that is different from the pressure profile of choles- 
terol systems. Unfortunately, no such simulations, let 
alone experimental measurements, are available at the 
moment. 



C. Centre of the bilayer 

The centre of the bilayer is formed by the ends of the 
acyl chains of the DPPC and cholesterol molecules. Both 
mass and electron density have a minimum at z = 
known as methyl trough. For low cholesterol concen- 
tration, there is a single peak in the pressure profile 
at z = that has been suggested to be related to in- 
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FIG. 4 Mass density profiles, with the contributions from 
both leaflets plotted separately (left for DPPC, right for 
cholesterol). The curves have been scaled by their maximum 
and shifted vertically. For symmetry reasons, only a small 
part of the right leaflet is shown. The numbers on the curves 
give the cholesterol concentration. The amount of interdigi- 
tation is significantly reduced if cholesterol is present. 



terdigitation ijvan den Brink- van der Laan et all [2004). 
Understanding the centre of the bilayer is difficult 
as the peak in the pressur e profile is found only in 
atom i stic MP simulations (iGullingsrud and Schulteii 
12004 iLindahl and Edholml l2000(l whereas analytical 
the ories and coarse-grained models f ail to reproduce 
it ijvan den Brink- van der Laan et"al 1 I2004D . Experi- 
mentally it was found that cholesterol reduces interdigita- 
tion but it should be noted that even pure DPPC bilayers 
are o nly moderately interdigitated ijSiminovitch et all 
This is in agreement with our results in Fig. 

For increasing cholesterol concentration, the single 
peak in the pressure profile at z — splits into two peaks 
as the pressure close to z = increases faster than pre- 
cisely at z = 0. Simultaneously, in addition to a reduc- 
tion in interdigitation, the end point of the acyl chains 
becomes more strongly defined pronounced in the mass 
density profile. At the highest cholesterol concentrations, 
even a small peak in the mass density profiles in Fig. 0] 
develops near the end of the chains. 

The relation between interdigitation and pressure pro- 
file in the centre thus is not obvious - our results are more 
consistent with the absence of such a relation. From 
Fig. |2| it is seen that the peak in the pressure profile 
around the bilayer centre does not originate in an in- 
crease of the expanding pressure components but rather 
from a sharp decrease of the condensing pressure. A pos- 
sible explanation is that, as they are less ordered than 
the other atoms, the atoms at the end of acyl chains can- 
not make favourable contact with neighbouring chains, 
thereby loosing van der Waals-attraction. The shape of 



the pressure profile near the centre of the bilayer would 
thus be dominated not by the mass density profile it- 
self but by the spatial distribution of end points of the 
chains. For cholesterol, this gives a peak at z = for 
small cholesterol concentration that splits into two peaks 
as interdigitation is reduced, in agreement with the com- 
puted pressure profiles. 



VI. RELATION TO MEMBRANE PENETRATION 

The lateral pressure profile is related to the abil- 
ity of small molecules to penetrate into the bilayer. 
Negative pressure means that the system would like 
to contract. Regions with negative lateral pressure 
thus mark regions where it would be energetically 
favourable for additional particles or molecules to be in- 
serted. Ultimately, the partitioning of a particular solute 
into the bilayer is determined by the free-energy pro- 
file ijMarrink and Berendsenl Il996j) . The lateral pressure 
profile gives the volume contribution to the free-energy 
profile, and for nonpolar or moderately polar solutes, this 
is a good approximation to the full free-energy profile. 

On first view, this seems to be related to free-volume 
theories which state that the amount of unoccupied space 
inside the bilayer determines its penetrability. In the con- 
text of pressure profiles, "occupied volume" is equivalent 
to large steric repulsive forces, resulting in a large ex- 
panding static pressure. Already from the discussion in 
Sec.[V]it should be obvious that free- volume theories face 
a big problem as much of the pressure inside the bilayer 
is either attractive or of kinetic origin. 

One of the many reasons why cholesterol contain- 
ing membranes are interesting systems to study is that 
cholester ol affects the membran e spatially inhomoge- 
neously ijLee and Petersenl |2004|) . For example, the ad- 
dition of cholesterol increases the penetration of water 
into most regions of the bilayer, while it reduces the wa- 
ter co ncentration in the very centre of the bilayer l|Marshl 
120021) . 

These observations are easily explained within the con- 
text of lateral pressure profiles but not within free- volume 
theories. In the pure DPPC bilayer, the pressure in the 
acyl chain region is small but positive. As the cholesterol 
concentration is increased from % to 30 %, the pressure 
in most of the acyl chain region becomes negative (al- 
lowing for easier penetration) while in the centre of the 
bilayer it becomes more positive (thus expelling particles 
there). Also at 50% cholesterol, there are ample regions 
with negative pressure even though they become discon- 
nected by high peaks of positive pressure. 

To make contact with free- volume theories, it is in- 
structive to "compute" from Fig. [21a "fake" pressure pro- 
file that only includes static expanding pressure contri- 
butions, shown in Fig. [BJ The "correct" lateral-pressure 
profile changes strongly as cholesterol is added, with sev- 
eral new peaks appearing in the pressure profile, and the 
magnitude of the pressure changing by about one order 
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FIG. 5 Pressure profiles for 30 % cholesterol. The lower curve 
gives the correct pressure profile whereas in the upper curve 
only static expanding contributions to the pressure have been 
included. 



of magnitude. In contrast, the free- volume profile lacks 
the additional structure induced by cholesterol and only 
changes monotonously as the cholesterol concentration is 
changed, and the changes are rather moder ate (at least 
compare d to the lateral-pressure profile) ijFalck et all 
l2004albt) - both observations apply also to the "fake" 
pressure profile. 

The free-volume fraction is highest in the centre of 
the bilayer, which is simply a refl ection of the bilaye r 
trough with reduced atom density ijFalck et all |2004a) , 
and molecules such as water should thus preferably be 
located there. From this contradiction, it was thus re- 
alised already early that it is not the average amount of 
free volume that determines the penetrat i on bu t rather 
rare fluctuations ijMarrink and Berendsenl Il994j) . 

Free-volume theory can be applied to compare pene- 
tration of solutes of different size through the same bi- 
layer, as penetration then is rate-limited by a single point 
in th e membrane, usually l ocate d at the bilayer inter- 
face llKillian and van Meerl 1200 1|) . Already when com- 
paring the same bilayer at different temper atures, such a 
one-p arameter approach may break down ijSutter et all 
2004). The lateral -pressure profile, on the other hand, 
allows to compare also very different systems, such as 
bilayers with different cholesterol concentration, as the 
description of the bilayer is in full detail. (The descrip- 
tion of the solute is still reduced to solely its volume, and 
the most straight-forward improvement to this simplifi- 
cation are free-energy profiles.) 



VII. DISCUSSION 

In this paper we have studied lateral-pressure profiles 
in cholesteroh-DPPC bilayers at varying composition ra- 



tios. We found that the pressure profiles become more 
structured and complicated as the cholesterol fraction 
is increased (cf. Fig. This is related to the struc- 
tural changes in the bilayer upon addition of cholesterol. 
These changes can have biological consequences: pro- 
teins or other molecules are able to anchor at a given 
dep th inside a cell membrane guided by the pr essure pro- 
file l)van den Brink- van der Laan et alll2004[) . 

Biological interactions can be divided into specific and 
unspecific. The general view in the literature seems to 
be that changes via the lateral-pressure profile are rather 
unspecific, i. e., proteins or other molecules in a cell mem- 
brane are affected by the lateral-pressure profile but it 
is of secondary importance of how a change in pressure- 
profile was induced. For specific interaction, a direct con- 
tact between cholesterol and, e.g., a membrane channel, 
would be needed. 

In view of the lateral-pressure profiles presented in this 
paper one could speculate whether a specific interaction 
via the pressure-profile is possible. The profile contains a 
wealth of additional structure, very different from what 
could be achieved by simply compressing the bilayer. 
This structure in the lateral-pressure profile might be 
specific enough to allo w molecules such as the nicotinic 
acetylcholine receptor (jRankin et all Il997(l to sense the 
presence or absence of cholesterol solely via the lateral- 
pressure profile 

This might be an additional reason for abundance of 
cholesterol in eukaryotic cell membranes, and why there 
is so much cholesterol but so few other chemically closely 
related steroids - in contrast to the multitude of different 
lipids found in cell membranes. From that point of view, 
it would be interesting to see the difference between the 
different sterols on pressure profiles, e.g., how the pres- 
sure profile would change if cholesterol was replaced by 
lanosterol. 

The results presented in this paper also demonstrate 
that the effects of cholesterol cannot be captured by sim- 
ple models. Still, generic models have their advantages, 
and the importance of lateral pressure profiles would 
not have been accepted without the p ioneering analyt- 
ical wor k by Cantor in the late 1990's ijCantort Il997albl 
Il999albf ). His most important contribution properly were 
not so much the results themselves - much of it was al- 
ready known - but to cast them into a single consis- 
tent form that, e. g., allowed to relate pressure profiles to 
"nonbilayer lipids" or spontaneous membrane curvature. 

Characteristic of such analytical theories is 
that they are very successful in describing plain 
acyl chains but already their treatment o f the 
headgroups is rather crude |Ben-Shaull Il995l 
Ivan den Brink- van der Laan et all l2004|) . Choles- 
terol is well beyond the reach of such models. In the 
headgroup region, the problem is the description of 
the mutual headgroup interactions and, even though 
this effect is generic, their change upon decrease in the 
area per lipid induced by cholesterol. In the acyl chain 
region, the problem is the non-generic nature of the rigid 
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cholesterol body. 

This and other atomistic studies show the power of 
computer simulations to resolve the details of molecu- 
lar systems and to provide further insight to complex 
multicomponent systems. Pressure profiles are a superb 
example for this power due to the severe limitations en- 
countered in experiment, coarse-grained simulations and 
analytical theories. 
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